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SOLUTION STRU CTURE OF TNFR-1 DP AND USES THEREOF 

5 This application claims the benefit of U.S. Provisional Application No. 

60/206,215 filed May 22, 2000. 

Field of the Invention 
The present invention relates to the three dimensional solution 
10 structure of tumor necrosis factor receptor 1 death domain (TNFR-1 DD), as well as 
the identification of various active site structures of TNFR-1 DD. These structures 
are critical for the design and selection of potent and selective inhibitors of TNF 
□ signaling pathways. 

: « 15 Background of the Invention 

^ Tumor Necrosis Factor Receptor 1 death domain (TNFR-1 DD) is the 

CP intracellular functional domain responsible for the receptor signaling activities. 

P Tumor necrosis factor receptor- 1, through trimerization induced by the binding of 

TNFa (cachectin) or TNFp (lymphotoxin a) trimers, elicits a variety of biological 
;F 20 responses including antiviral activity, cytotoxicity, and modulation of gene 
s £ expression [1, 2]. TNFa and lymphotoxin a also bind and activate TNFR-2 [3]. 

Based on similarities in their cysteine-rich extracellular domains TNFR-1 and TNFR- 
2 belong to a receptor superfamily, which besides a number of death inducing 
receptors, includes CD40 and the low-affinity nerve growth factor receptor [4] . 
25 Although most cell types express both TNF receptors, TNFR-1 appears to play a 
predominant role in induction of gene expression and induction of cell death by 
TNFa [5, 6]. 

The death domain was originally described as a region of similarity 
within the intracellular portions of the TNFR-1 and Fas/Apol that is essential for 
30 the transduction of cytotoxic signals [7] . It was subsequently shown that the death 
domain is a protein interaction motif involved in homo- and hetero-association [8] . 
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TNFR-1 trimerizes upon binding to the TNF trimer inducing association of its 
intracellular death domain [4] . The trimerization of TNFR-1 allows for the 
recruitment of an adaptor protein named TNFR-Associated Factor-2 (TRAF-2) [11] 
through interactions with the N-terminal domain of Fas Associated Death Domain 

5 (FADD) [11] through death domain interactions. 

The solution structure of Fas, FADD and p75 neurotrophin receptor 
death domains have recently been determined by NMR spectroscopy [12-14]. The 
overall fold observed for these death domains consist of six antiparallel 
amphipathic cc-helices packed into a globular structure. Two other domains called 

10 DED (Death effector domain) and CARD (Caspase recruitment domain) are 

structurally related to the death domain [15]. Recently the solution structure of 
the FADD-DED and RAIDD-CARD were solved revealing topologies very similar to 
the known death domain structures [16, 17]. Two recent studies have described 
the three-dimensional structure of Apaf-1 CARD in complex with the procaspase-9 

15 CARD [18] and Tube death domain in complex with Pelle death domain [19] . 
These different studies illustrate the diverse association mechanisms in the death 

domain superfamily. 

The low level of sequence conservation between death domains 
probably reflects their role in diverse cellular functions. Additionally, the apparent 

20 absence of a conserved interaction surface suggests that death domains may 
associate by a variety of mechanisms. Indeed, the surface formed by a-helices 2 
and 3 has been implicated in the homo- and hetero-association of the death 
domains of FAS [12] and FADD [13]. Similar observations have been made in the 
case of the interaction of the CARD domains of Apaf-1 and procaspase-9 [18, 20- 

25 22] . Conversely, the complex formation of the Tube and Pelle death domains 

seems to rely on contacts between a-helices 4 and 5 of the Pelle death domain with 
cc-helix 6 and the unique C-terminal tail of the Tube death domain [19] . The 
nature of the interaction also seems to vary between death domain complexes. 
Electrostatic interactions are thought to be a key component in the interaction 
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between FAS death domain (FAS DD) and FADD death domain (FADD DD) [12]. 
Whereas, hydrogen bond contacts and van der Waals interactions have been shown 
to be involved in the complex of Tube and Pelle death domains [19], as well as in 
the complex of Apaf-1 and procaspase-9 CARD domains [18]. It is currently 
believed that the interactions between DEDs are hydrophobic in nature [16]. 

Summary of the Invention 

The present invention relates to the three dimensional structure of a 
tumor necrosis factor receptor 1 death domain (TNFR-1 DD), and more specifically, 
to the solution structure of TNFR-1 DD, as determined using spectroscopy and 
various computer modeling techniques. 

Particularly, the invention is further directed to the identification, 
characterization and three dimensional structure of an active site of TNFR-1 DD 
that provides an attractive target for the rational design of potent and selective 
inhibitors of TNF signaling pathways. 

Accordingly, the present invention provides a solution comprising a 
tumor necrosis factor receptor 1 death domain (TNFR-1 DD). The three 
dimensional solution structure of TNFR-1 DD is provided by the relative atomic 
structural coordinates of Figure 8, as obtained from spectroscopy data. 

Also provided by the present invention is an active site of a TRADD DD 
binding protein or peptide, preferably of TNFR-1 DD, wherein said active site is 
characterized by a three dimensional structure comprising the relative structural 
coordinates of amino acid residues K343, E344, R347, R348 and D353 according to 
Figure 8, ± a root mean square deviation from the conserved backbone atoms of 

o 

said amino acids of not more than 1.5 A. 

Also provided for by the present invention is an active site of a TRADD 
DD binding protein or peptide, preferably of TNFR-1 DD, wherein said active site is 
characterized by a three dimensional structure comprising the relative structural 
coordinates of amino acid residues E369 and Y373 according to Figure 8, ± a root 
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mean square deviation from the conserved backbone atoms of said amino acids of 

o 

not more than 1.5 A. 

In addition, the present invention provides an active site of a TRADD 
DD binding protein or peptide, preferably TNFR-1 DD, wherein said active site is 
5 characterized by a three dimensional structure comprising the relative structural 
coordinates of amino acid residues E335, E386, E390, D398, E406, D407, E409 
and E410 according to Figure 8, ± a root mean square deviation from the 

o 

conserved backbone atoms of said amino acids of not more than 1.5 A. 

Still further, the present invention provides an active site of a TRADD 
10 DD binding protein or peptide, preferably TNFR-1 DD, wherein said active site is 
characterized by a three dimensional structure comprising the relative structural 
coordinates of amino acid residues R358, R365, R368, R379, R380, R381, R384, 
R385, R394 and R397 according to Figure 8, ± a root mean square deviation from 

o 

the conserved backbone atoms of said amino acids of not more than 1.5 A. 

15 The solution coordinates of TNFR-1 DD or portions thereof (such as the 

TRADD DD binding site or other active sites), as provided by this invention may be 
stored in a machine-readable form on a machine-readable storage medium, e.g. a 
computer hard drive, diskette, DAT tape, etc., for display as a three-dimensional 
shape or for other uses involving computer-assisted manipulation of, or 

20 computation based on, the structural coordinates or the three-dimensional 

structures they define. By way of example, the data defining the three dimensional 
structure of TNFR-1 DD as set forth in Figure 8 may be stored in a 
machine-readable storage medium, and may be displayed as a graphical 
three-dimensional representation of the relevant structural coordinates, typically 

25 using a computer capable of reading the data from said storage medium and 
programmed with instructions for creating the representation from such data. 

Accordingly, the present invention provides a machine, such as a 
computer, programmed in memory with the coordinates of TNFR-1 DD or portions 
thereof, together with a program capable of converting the coordinates into a three 
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dimensional graphical representation of the structural coordinates on a display 
connected to the machine. A machine having a memory containing such data aids 
in the rational design or selection of inhibitors of TNFR-1 DD binding or activity, 
including the evaluation of the ability of a particular chemical entity to favorably 
5 associate with TNFR-1 DD as disclosed herein, as well as in the modeling of 

compounds, proteins, complexes, etc. related by structural or sequence homology 
to TNFR-1 DD. 

The present invention is additionally directed to a method of 
determining the three dimensional structure of a molecule or molecular complex 

10 whose structure is unknown, comprising the steps of first obtaining crystals or a 
solution of the molecule or molecular complex whose structure is unknown, and 
then generating X-ray diffraction data from the crystallized molecule or molecular 
complex and/or generating NMR data from the solution of the molecule or 
molecular complex. The generated diffraction or spectroscopy data from the 

15 molecule or molecular complex can then be compared with the solution 

coordinates or three dimensional structure of TNFR-1 DD as disclosed herein, and 
the three dimensional structure of the unknown molecule or molecular complex 
conformed to the TNFR-1 DD structure using standard techniques such as 
molecular replacement analysis, 2D, 3D and 4D isotope filtering, editing and triple 

20 resonance NMR techniques, and computer homology modeling. Alternatively, a 
three dimensional model of the unknown molecule may be generated by generating 
a sequence alignment between TNFR-1 DD and the unknown molecule, based on 
any or all of amino acid sequence identity, secondary structure elements or tertiary 
folds, and then generating by computer modeling a three dimensional structure for 

25 the molecule using the three dimensional structure of, and sequence alignment 

with, TNFR-1 DD. 

The present invention further provides a method for identifying an 
agent that interacts with TNFR-1 DD, comprising the steps of determining an active 
site of TNFR-1 DD using the three dimensional TNFR-1 DD structure, and then 
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performing computer fitting analyses to identify an agent which interacts with the 
identified active site. Once the agent has been identified, it may be contacted with 
TNFR-1 DD (using TNFR-1 DD or a molecule comprising TNFR-1 DD such as TNFR- 
1), and the effect the agent has on TNFR-1 DD may then be assessed. In addition, 
5 the agent may be contacted with TNFR-1 DD (using TNFR-1 DD or a molecule 
comprising TNFR-1 DD such as TNFR-1) in the presence of a TNFR-1 DD binding 
molecule (including but not limited to TRADD DD), and the effect the agent has on 
binding between TNFR-1 DD and the TNFR-1 DD binding molecule may then be 
assessed. 

10 Also provided is a method for identifying a potential inhibitor of TNFR- 

1 DD, comprising the steps of using a three dimensional structure of TNFR-1 DD as 
defined by the relative structural coordinates of amino acids encoding TNFR-1 DD 
to design or select a potential inhibitor, and obtaining or synthesizing said potential 
inhibitor. The inhibitor may be selected by screening an appropriate database, may 

15 designed de novo by analyzing the steric configurations and charge potentials of an 
empty TNFR-1 DD active site in conjunction with the appropriate software 
programs, or may be designed using characteristics of known inhibitors to create 
"hybrid" inhibitors. The inhibitor may then be contacted with TNFR-1 DD alone 
(using TNFR-1 DD or a molecule comprising TNFR-1 DD such as TNFR-1), or in the 

20 presence of a TNFR-1 DD binding molecule such as TRADD DD, and the effect of 
the inhibitor on TNFR-1 alone or binding between TNFR-1 and the TNFR-1 DD 
binding molecule may be assessed. It is also within the confines of the present 
invention that a potential inhibitor may be designed or selected by identifying 
chemical entities or fragments capable of associating with TNFR-1 DD; and 

25 assembling the identified chemical entities or fragments into a single molecule to 
provide the structure of the potential inhibitor. 

Finally, the present invention provides agents or inhibitors designed or 
selected using the methods disclosed herein. Additional objects of the present 
invention will be apparent from the description which follows. 
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Brief Description of the Figures 
Figure 1 is a summary of NMR-derived secondary structure indicators 
measured for TNFR-DD R347A together with the deduced secondary structure. For 
5 NOE data, the thickness of the lines reflects strength of sequential NOEs. 

Figure 2 is a summary of the inter-helical NOE contacts of TNFR-DD 
R347A at pH8.8. The helices 1, 3 and 5 are viewed from their N-termini, while 
helices 2, 4 and 6 are viewed from their C-termini. The bold and boxed residues 
have altered binding properties when mutated to alanine (see text and Fig. 5) . The 
10 bold and italic residues have unaltered binding properties when mutated to alanine 
(see text and Fig. 5) . 

Figure 3 depicts the effect of the ionic strength on the binding of TNFR- 
DD (200ng/ml) to MBP-TNFR-DD (left bar), MBP-TRADD (right bar). 

Figure 4 represents a sequence alignment of the death domains of 
15 TNFR-1 and Fas. The helices in the NMR structure of the TNFR-1 and Fas (Huang 
et al., 1996) are indicated. The mutations are indicated by an arrow pointing to 
the amino acid of substitution. A minus sign above the TNFR-1 mutation indicates 
a loss in self-association or interaction with TRADD-DD (see text). A minus sign 
below the Fas mutation indicates that this mutant aggregates less than the wild- 
20 type protein (Huang et al, 1996). 

Figures 5A and 5B depict an ELISA of TNFR-DD wild-type and mutants 
binding to MBP-TNFR-DD and MBP-TRADD. Figure 5A: Binding of the TNFR-DD 
mutants to MBP-TNFR-DD immobilized on the plate at 200 ng/well. The reference 
at 100% is the binding of the wild-type protein to MBP-TNFR-DD. The values used 
25 for the determination of the binding percentage were obtained with -300 ng/ml of 
soluble protein. Figure 5B: Binding of the TNFR-DD mutants to MBP-TRADD. The 
experimental conditions used were the same as with MBP-TNFR-DD. 

Figure 6 depicts a summary of H-D exchange, 3 J HNHK scalar coupling 
information, sequential NOEs (d aN , d NN , dp N ), medium range NOEs (d NN(i i+2) , 
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<Wi +3 ))> d aNai+4> d aPai+3] ) 3 and the A 13 C a and A 13 C p . Amides that didn't exchange 

with D 2 0 within 30 minutes are designated with an (•). The residues having a 
3 J HNHC£ less than 5 Hz are designated with an solid line. The intensities of the NOEs 
are represented by the thickness of the lines. The values of the A 13 C a and A 13 C p are 
5 represented by the intensity of the blocks. The secondary structure highlighting the 
6 a-helices is indicated at the bottom. 

Figure 7 depicts the secondary structure and sequence alignment of 
TNFR-1 DD with proteins in the death domain superfamily. In instances when the 
lengths of the helices differed when compared to TNFR-1 DD, the sequences were 
10 aligned based upon sequence homology instead of secondary structure alignment. 

Figure 8 lists the atomic structure coordinates for the restrained 
minimized mean structure of TNFR-1 DD as derived by multidimensional NMR 
spectroscopy. "Atom type" refers to the atom whose coordinates are being 
measured. "Residue" refers to the type of residue of which each measured atom is a 
15 part - i.e., amino acid, cofactor, ligand or solvent. The "x, y and z" coordinates 

o 

indicate the Cartesian coordinates of each measured atom's location (A) . All non- 
protein atoms are listed as HETATM instead of atoms using PDB conventions. 



Detailed Description of the Invention 
20 As used herein, the following terms and phrases shall have the 

meanings set forth below: 

Unless otherwise noted, "TNFR-1 DD" includes both the death domain 
of TNFR-1 as encoded by the amino acid sequence of Figure 6 (including 
conservative substitutions thereof), as well as "TNFR-1 DD analogues", defined 
25 herein as proteins or peptides comprising a TRADD DD or TRADD DD-like binding 
site. Such TNFR-1 DD analogues include, but are not limited to, an active site 
characterized by a three dimensional structure comprising (i) the relative structural 
coordinates of amino acid residues K343, E344, R347, R348 and D353 according to 
Figure 8, (ii) the relative structural coordinates of amino acid residues E369 and 
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Y373 according to Figure 8, (iii) the relative structural coordinates of amino acid 
residues E335, E386, E390, D398, E406, D407, E409 and E410 according to Figure 
8, or (iv) the relative structural coordinates of amino acid residues R358, R365, 
R368, R379, R380, R381, R384, R385, R394 and R397 according to Figure 8, in 
5 each case, ± a root mean square deviation from the conserved backbone atoms of 

o o 

said amino acids of not more than 1.5 A, more preferably not more than 1.0 A, and 

o 

most preferably not more than 0.5 A. 

Unless otherwise indicated, "protein" or "molecule" shall include a 
protein, protein domain, polypeptide or peptide. 

10 "Structural coordinates" are the Cartesian coordinates corresponding to 

an atom's spatial relationship to other atoms in a molecule or molecular complex. 
Structural coordinates may be obtained using x-ray crystallography techniques or 
NMR techniques, or may be derived using molecular replacement analysis or 
homology modeling. Various software programs allow for the graphical 

15 representation of a set of structural coordinates to obtain a three dimensional 

representation of a molecule or molecular complex. The structural coordinates of 
the present invention may be modified from the original set provided in Figure 8 by 
mathematical manipulation, such as by inversion or integer additions or 
subtractions. As such, it is recognized that the structural coordinates of the present 

20 invention are relative, and are in no way specifically limited by the actual x, y, z 
coordinates of Figure 8. 

An "agent" shall include a protein, polypeptide, peptide, nucleic acid, 
including DNA or RNA, molecule, compound or drug. 

"Root mean square deviation" is the square root of the arithmetic mean 

25 of the squares of the deviations from the mean, and is a way of expressing 
deviation or variation from the structural coordinates described herein. The 
present invention includes all embodiments comprising conservative substitutions 
of the noted amino acid residues resulting in same structural coordinates within the 
stated root mean square deviation. 
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It will be obvious to the skilled practitioner that the numbering of the 
amino acid residues in the various isoforms of TNFR-1 DD covered by the present 
invention may be different than that set forth herein, or may contain certain 
conservative amino acid substitutions that yield the same three dimensional or 
5 solution structures as those defined by Figure 8 herein. Corresponding amino acids 
and conservative substitutions in other isoforms or analogues are easily identified 
by visual inspection of the relevant amino acid sequences or by using commercially 
available homology software programs. 

"Conservative substitutions" are those amino acid substitutions which 

10 are functionally equivalent to the substituted amino acid residue, either by way of 
having similar polarity, steric arrangement, or by belonging to the same class as the 
substituted residue (e.g., hydrophobic, acidic or basic), and includes substitutions 
having an inconsequential effect on the three dimensional structure of TNFR-1 DD 
with respect to the use of said structure for the identification and design of TNFR-1 

15 DD or TNFR-1 DD complex inhibitors, for molecular replacement analyses and/or 
for homology modeling. 

An "active site" refers to a region of a molecule or molecular complex 
that, as a result of its shape and charge potential, favorably interacts or associates 
with another agent (including, without limitation, a protein, polypeptide, peptide, 

20 nucleic acid, including DNA or RNA, molecule, compound, antibiotic or drug) via 
various covalent and/or non-covalent binding forces. As such, an active site of the 
present invention may include, for example, the actual site of TRADD DD binding 
with TNFR-1 DD, as well as accessory binding sites adjacent to the actual site of 
TRADD DD binding that nonetheless may affect TNFR-1 DD upon interaction or 

25 association with a particular agent, either by direct interference with the actual site 
of TRADD DD binding or by indirectly affecting the steric conformation or charge 
potential of TNFR-1 DD and thereby preventing or reducing TRADD DD binding to 
TNFR-1 DD at the actual site of TRADD DD binding. As used herein, "active site" 
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also includes the TNFR-1 DD site of self association, as well as other binding sites 
present on TNFR-1 DD. 

A "TNFR-1 DD complex" refers to a co-complex of a molecule 
comprising the TNFR-1 DD region in bound association with a protein, polypeptide, 
5 peptide, nucleic acid, including DNA or RNA, small molecule, compound or drug, 
either by covalent or non-covalent binding forces. A non-limiting example of a 
TNFR-1 DD complex includes TNFR-1, TNFR-1 DD or a TNFR-1 DD analogue 
bound to TRADD DD. 

The present invention relates to the three dimensional structure of 
10 TNFR-1 DD or of a TNFR-1 DD analogue, and more specifically, to the solution 
O structure of TNFR-1 DD as determined using multidimensional NMR spectroscopy 

i;0 and various computer modeling techniques. The solution coordinates of TNFR-1 

DD (disclosed herein at Figure 8) are useful for a number of applications, 
;5 including, but not limited to, the characterization of a three dimensional structure 

W 15 of TNFR-1 DD, as well as the visualization, identification and characterization of 

□ TNFR-1 DD active sites, including the site of TRADD DD binding to TNFR-1 DD. 

!^ The active site structures may then be used to predict the orientation and binding 

J affinity of a designed or selected inhibitor of TNFR-1 DD, a TNFR-1 DD analogue or 

of a TNFR-1 DD complex. Accordingly, the invention is particularly directed to the 
20 three dimensional structure of a TNFR-1 DD active site, including but not limited to 
the TRADD DD binding site. 

As used herein, the TNFR-1 DD in solution comprises amino acid 
residues 316-425 of TNFR-1, and more preferably, the amino acid residues 316- 
425 set forth in Figure 6, or conservative substitutions thereof. Preferably, the 
25 TNFR-1 DD in solution is either unlabeled, 15 N enriched or 15 N, 13 C enriched, and is 
preferably biologically active. In addition, the secondary structure of the TNFR-1 
DD in the solutions of the present invention comprises six alpha (a) helices. In this 
regard, al comprises amino acid residues A328-N336 of TNFR-1 DD, a2 comprises 
amino acid residues W342-L349 of TNFR-1 DD, 063 comprises amino acid residues 
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P353-L361 of TNFR-1 DD, a4 comprises amino acid residues L367-R380 of TNFR-1 
DD, a5 comprises amino acid residues L389-D398 of TNFR-1 DD and a6 comprises 
amino acid residues G403-L412 of TNFR-1 DD. 

The protein or peptide used in the solution of the present invention 
5 includes TNFR-1 DD, as well as TNFR-1 DD analogues, where said protein or 
peptide preferably comprises an active site characterized by a three dimensional 
structure comprising (i) the relative structural coordinates of amino acid residues 
K343, E344, R347, R348 and D353 according to Figure 8, (ii) the relative 
structural coordinates of amino acid residues E369 and Y373 according to Figure 8, 
10 (iii) the relative structural coordinates of amino acid residues E335, E386, E390, 
D398, E406, D407, E409 and E410 according to Figure 8, or (iv) the relative 
structural coordinates of amino acid residues R358, R365, R368, R379, R380, 
R381, R384, R385, R394 and R397 according to Figure 8, in each case, ± a root 
mean square deviation from the conserved backbone atoms of said amino acids of 

o o 

15 not more than 1.5 A, more preferably not more than 1.0 A, and most preferably not 

o 

more than 0.5A. In the most preferred embodiment, the protein or peptide used in 
the solution of the present invention is characterized by a three dimensional 
structure comprising the complete structural coordinates of the amino acids 
according to Figure 8, ± a root mean square deviation from the conserved 

o 

20 backbone atoms of said amino acids of not more than 1.5 A (or more preferably, 

o o 

not more than 1.0 A, and most preferably, not more than 0.5 A). 

Molecular modeling methods known in the art may be used to identify 
an active site or binding pocket of TNFR-1 DD, a TNFR-1 DD complex, or of a 
TNFR-1 DD analogue. Specifically, the solution structural coordinates provided by 
25 the present invention may be used to characterize a three dimensional structure of 
the TNFR-1 DD molecule, molecular complex or TNFR-1 DD analogue. From such 
a structure, putative active sites may be computationally visualized, identified and 
characterized based on the surface structure of the molecule, surface charge, steric 
arrangement, the presence of reactive amino acids, regions of hydrophobicity or 
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hydrophilicity, etc. Such putative active sites may be further refined using chemical 
shift perturbations of spectra generated from various and distinct TNFR-1 DD 
complexes, competitive and non-competitive inhibition experiments, and/or by the 
generation and characterization of TNFR-1 DD or ligand mutants to identify critical 
5 residues or characteristics of the active site. 

The identification of putative active sites of a molecule or molecular 
complex is of great importance, as most often the biological activity of a molecule 
or molecular complex results from the interaction between an agent and one or 
more active sites of the molecule or molecular complex. Accordingly, the active 
10 sites of a molecule or molecular complex are the best targets to use in the design or 
O selection of inhibitors that affect the activity of the molecule or molecular complex. 

iM. The present invention is directed to an active site of TNFR 1 DD, a 

lp TNFR-1 DD complex or of a TNFR-1 DD analogue, that, as a result of its shape, 

reactivity, charge potential, etc., favorably interacts or associates with another 
isM 15 agent (including, without limitation, a protein, polypeptide, peptide, nucleic acid, 

□ including DNA or RNA, molecule, compound, antibiotic or drug) . Preferably, the 

u the present invention is directed to an active site of a TRADD DD binding protein or 

;:£ peptide, and preferably TNFR-1 DD, that is characterized by the three dimensional 

H structure comprising the relative structural coordinates of amino acid residues 

20 K343, E344, R347, R348 and D353 according to Figure 8, ± a root mean square 
deviation from the conserved backbone atoms of said amino acids of not more than 

o o o 

1.5 A, preferably not more than 1.0 A, and most preferably not more than 0.5 A. 
In another embodiment, the active site of a TRADD DD binding protein or peptide, 
and preferably TNFR-1 DD, is characterized by the three dimensional structure 
25 comprising the relative structural coordinates of amino acid residues E369 and 
Y373 according to Figure 8, ± a root mean square deviation from the conserved 

o 

backbone atoms of said amino acids of not more than 1.5 A, preferably not more 

o o 

than 1.0 A, and most preferably not more than 0.5 A. In another embodiment, the 
active site of a TRADD DD binding protein or peptide, and preferably TNFR-1 DD, 
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is characterized by the three dimensional structure comprising the relative 
structural coordinates of amino acid residues E335, E386, E390, D398, E406, 
D407, E409 and E410 according to Figure 8, ± a root mean square deviation from 

o 

the conserved backbone atoms of said amino acids of not more than 1.5 A, 

o o 

5 preferably not more than 1.0 A, and most preferably not more than 0.5 A. In yet 
another embodiment, the active site of a TRADD DD binding protein or peptide, 
and preferably TNFR-1 DD, is characterized by the three dimensional structure 
comprising the relative structural coordinates of amino acid residues R358, R365, 
R368, R379, R380, R381, R384, R385, R394 and R397 according to Figure 8, ± a 
10 root mean square deviation from the conserved backbone atoms of said amino 

o o 

acids of not more than 1.5 A, preferably not more than 1.0 A, and most preferably 
not more than 0.5 A. 

In order to use the structural coordinates generated for a solution 
structure of the present invention as set forth in Figure 8, it is often necessary to 

15 display the relevant coordinates as, or convert them to, a three dimensional shape 
or graphical representation, or to otherwise manipulate them. For example, a three 
dimensional representation of the structural coordinates is often used in rational 
drug design, molecular replacement analysis, homology modeling, and mutation 
analysis. This is typically accomplished using any of a wide variety of commercially 

20 available software programs capable of generating three dimensional graphical 
representations of molecules or portions thereof from a set of structural 
coordinates. Examples of said commercially available software programs include, 
without limitation, the following: GRID (Oxford University, Oxford, UK); MCSS 
(Molecular Simulations, San Diego, CA); AUTODOCK (Scripps Research Institute, 

25 La Jolla, CA); DOCK (University of California, San Francisco, CA); Flo99 

(Thistlesoft, Morris Township, NJ); Ludi (Molecular Simulations, San Diego, CA); 
QUANTA (Molecular Simulations, San Diego, CA); Insight (Molecular Simulations, 
San Diego, CA); SYBYL (TRIPOS, Inc., St. Louis. MO); and LEAPFROG (TRIPOS, 
Inc., St. Louis, MO). 
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For storage, transfer and use with such programs, a machine, such as a 
computer, is provided for that produces a three dimensional representation of the 
TNFR-1 DD, a portion thereof (such as an active site or a binding site), a TNFR-1 
DD molecular complex, or a TNFR-1 DD analogue. The machine of the present 
5 invention comprises a machine-readable data storage medium comprising a data 
storage material encoded with machine-readable data. Machine-readable storage 
media comprising data storage material include conventional computer hard drives, 
floppy disks, DAT tape, CD-ROM, and other magnetic, magneto-optical, optical, 
floptical and other media which may be adapted for use with a computer. The 

10 machine of the present invention also comprises a working memory for storing 
instructions for processing the machine-readable data, as well as a central 
processing unit (CPU) coupled to the working memory and to the machine- 
readable data storage medium for the purpose of processing the machine-readable 
data into the desired three dimensional representation. Finally, the machine of the 

15 present invention further comprises a display connected to the CPU so that the 

three dimensional representation may be visualized by the user. Accordingly, when 
used with a machine programmed with instructions for using said data, e.g., a 
computer loaded with one or more programs of the sort identified above, the 
machine provided for herein is capable of displaying a graphical three-dimensional 

20 representation of any of the molecules or molecular complexes, or portions of 
molecules of molecular complexes, described herein. 

In one embodiment of the invention, the machine-readable data 
comprises the relative structural coordinates of amino acid residues K343, E344, 
R347, R348 and D353, ± a root mean square deviation from the conserved 

o 

25 backbone atoms of said amino acids of not more than 1.5 A, or preferably, not 

o o 

more than 1.0 A, or more preferably not more than 0.5 A. In an alternate 
embodiment, the machine-readable data further comprises the relative structural 
coordinates of amino acid residues E369 and Y373 according to Figure 8, ± a root 
mean square deviation from the conserved backbone atoms of said amino acids of 
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not more than 1.5 A, preferably not more than 1.0 A, and most preferably not 

0 

more than 0.5 A. In another embodiment, the machine-readable data still further 
comprises the relative structural coordinates of amino acid residues E335, E386, 
E390, D398, E406, D407, E409 and E410 according to Figure 8, ± a root mean 
5 square deviation from the conserved backbone atoms of said amino acids of not 

0 0 

more than 1.5 A, or preferably, not more than 1.0 A, or more preferably not more 

o 

than 0.5 A. In yet another embodiment, the machine readable data comprises the 
relative structural coordinates of amino acid residues R358, R365, R368, R379, 
R380, R381, R384, R385, R394 and R397 according to Figure 8, ± a root mean 
10 square deviation from the conserved backbone atoms of said amino acids of not 

o o 

more than 1.5 A, or preferably, not more than 1.0 A, or more preferably not more 

0 

than 0.5 A. Finally, in the most preferred embodiment, the machine readable data 
comprises the complete structural coordinates according to Figure 8, ± a root mean 

o o 

square deviation of not more than 1.5 A (or more preferably, not more than 1.0 A, 

o 

15 and most preferably, not more than 0.5 A). 

The structural coordinates of the present invention permit the use of 
various molecular design and analysis techniques in order to (i) solve the three 
dimensional structures of related molecules, molecular complexes or TNFR-1 DD 
analogues, and (ii) to design, select, and synthesize chemical agents capable of 

20 favorably associating or interacting with an active site of an TNFR-1 DD molecule, 
molecular complex or TNFR-1 DD analogue, wherein said chemical agents 
potentially act as inhibitors of TNFR-1 DD or TNFR-1 DD complex binding to a 
number of binding proteins, including, but not limited to, TRADD DD. 

More specifically, the present invention provides a method for 

25 determining the molecular structure of a molecule or molecular complex whose 
structure is unknown, comprising the steps of obtaining a solution of the molecule 
or molecular complex whose structure is unknown, and then generating NMR data 
from the solution of the molecule or molecular complex. The NMR data from the 
molecule or molecular complex whose structure is unknown is then compared to 
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the solution structure data obtained from the TNFR-1 DD solutions of the present 
invention. Then, 2D, 3D and 4D isotope filtering, editing and triple resonance 
NMR techniques are used to conform the three dimensional structure determined 
from the TNFR-1 DD solution of the present invention to the NMR data from the 
5 solution molecule or molecular complex. Alternatively, molecular replacement may 
be used to conform the TNFR-1 DD solution structure of the present invention to x- 
ray diffraction data from crystals of the unknown molecule or molecular complex. 

Molecular replacement uses a molecule having a known structure as a 
starting point to model the structure of an unknown crystalline sample. This 
10 technique is based on the principle that two molecules which have similar 

structures, orientations and positions will diffract x-rays similarly. A corresponding 
approach to molecular replacement is applicable to modeling an unknown solution 
structure using NMR technology. The NMR spectra and resulting analysis of the 
NMR data for two similar structures will be essentially identical for regions of the 
15 proteins that are structurally conserved, where the NMR analysis consists of 
obtaining the NMR resonance assignments and the structural constraint 
assignments, which may contain hydrogen bond, distance, dihedral angle, coupling 
constant, chemical shift and dipolar coupling constant constraints. The observed 
differences in the NMR spectra of the two structures will highlight the differences 
20 between the two structures and identify the corresponding differences in the 
structural constraints. The structure determination process for the unknown 
structure is then based on modifying the NMR constraints from the known structure 
to be consistent with the observed spectral differences between the NMR spectra. 

Accordingly, in one non-limiting embodiment of the invention, the 
25 resonance assignments for the TNFR-1 DD solution provide the starting point for 
resonance assignments of TNFR-1 DD in a new TNFR-1 DD:"unsolved agent" 
complex. Chemical shift perturbances in two dimensional ls N/ l H spectra can be 
observed and compared between the TNFR-1 DD solution and the new TNFR-1 
DD: agent complex. In this way, the affected residues may be correlated with the 
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three dimensional structure of TNFR-1 DD as provided by the relevant structural 
coordinates of Figure 8. This effectively identifies the region of the TNFR-1 
DD: agent complex that has incurred a structural change relative to the native 
TNFR-1 DD structure. The a H, 15 N, 13 C and 13 CO NMR resonance assignments 
corresponding to both the sequential backbone and side-chain amino acid 
assignments of TNFR-1 DD may then be obtained and the three dimensional 
structure of the new TNFR-1 DD: agent complex may be generated using standard 
2D, 3D and 4D triple resonance NMR techniques and NMR assignment 
methodology, using the TNFR-1 DD solution structure, resonance assignments and 
structural constraints as a reference. Various computer fitting analyses of the new 
agent with the three dimensional model of TNFR-1 DD may be performed in order 
to generate an initial three dimensional model of the new agent complexed with 
TNFR-1 DD, and the resulting three dimensional model may be refined using 
standard experimental constraints and energy minimization techniques in order to 
position and orient the new agent in association with the three dimensional 

structure of TNFR-1 DD. 

The present invention further provides that the structural coordinates of 
the present invention may be used with standard homology modeling techniques 
in order to determine the unknown three-dimensional structure of a molecule or 
molecular complex. Homology modeling involves constructing a model of an 
unknown structure using structural coordinates of one or more related protein 
molecules, molecular complexes or parts thereof (i.e., active sites). Homology 
modeling may be conducted by fitting common or homologous portions of the 
protein whose three dimensional structure is to be solved to the three dimensional 
structure of homologous structural elements in the known molecule, specifically 
using the relevant (i.e., homologous) structural coordinates provided by Figure 8 
herein. Homology may be determined using amino acid sequence identity, 
homologous secondary structure elements, and/or homologous tertiary folds. 
Homology modeling can include rebuilding part or all of a three dimensional 
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structure with replacement of amino acids (or other components) by those of the 
related structure to be solved. 

Accordingly, a three dimensional structure for the unknown molecule 
or molecular complex may be generated using the three dimensional structure of 
5 the TNFR-1 DD molecule of the present invention, refined using a number of 
techniques well known in the art, and then used in the same fashion as the 
structural coordinates of the present invention, for instance, in applications 
involving molecular replacement analysis, homology modeling, and rational drug 
design. 

10 Determination of the three dimensional structure of TNFR-1 DD, its 

^ TRADD DD binding active site, and other binding sites, is critical to the rational 

W identification and/or design of agents that may act as inhibitors of TNFR-1 DD, 

i such as inhibitors of TRADD DD binding to TNFR-1 DD. This is advantageous over 

q conventional drug assay techniques, in which the only way to identify such an 

y 15 agent is to screen thousands of test compounds until an agent having the desired 

□ inhibitory effect on a target compound is identified. Necessarily, such conventional 

U screening methods are expensive, time consuming, and do not elucidate the 

S method of action of the identified agent on the target compound. 

^ : However, advancing X-ray, spectroscopic and computer modeling 

20 technologies allow researchers to visualize the three dimensional structure of a 
targeted compound (i.e., of TNFR-1 DD). Using such a three dimensional 
structure, researchers identify putative binding sites and then identify or design 
agents to interact with these binding sites. These agents are then screened for an 
inhibitory effect upon the target molecule. In this manner, not only are the number 
25 of agents to be screened for the desired activity greatly reduced, but the mechanism 
of action on the target compound is better understood. 

Accordingly, the present invention further provides a method for 
identifying a potential inhibitor of TNFR-1 DD, a TNFR-1 DD analogue or of a 
TNFR-1 DD complex, comprising the steps of using a three dimensional structure of 
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TNFR-1 DD as defined by the relative structural coordinates of Figure 8 to design or 
select a potential inhibitor of TNFR-1 DD activity, and synthesizing or obtaining 
said potential inhibitor. The inhibitor may be selected by screening an appropriate 
database, may be designed de novo by analyzing the steric configurations and 

5 charge potentials of an empty TNFR-1 DD or TNFR-1 DD complex active site in 
conjunction with the appropriate software programs, or may be designed using 
characteristics of known inhibitors of protein binding to TNFR-1 DD or TNFR-1 DD 
complexes in order to create "hybrid" inhibitors. 

An agent that interacts or associates with an active site of TNFR-1 DD, a 

10 TNFR-1 DD complex or a TNFR-1 DD analogue may be identified by determining 
an active site from the three dimensional structure of TNFR-1 DD, and performing 
computer fitting analyses to identify an agent which interacts or associates with 
said active site. Computer fitting analyses utilize various computer software 
programs that evaluate the "fit" between the putative active site and the identified 

15 agent, by (a) generating a three dimensional model of the putative active site of a 
molecule or molecular complex using homology modeling or the atomic structural 
coordinates of the active site, and (b) determining the degree of association 
between the putative active site and the identified agent. The degree of association 
may be determined computationally by any number of commercially available 

20 software programs, or may be determined experimentally using standard binding 
assays. 

Three dimensional models of the putative active site may be generated 
using any one of a number of methods known in the art, and include, but are not 
limited to, homology modeling as well as computer analysis of raw structural 
25 coordinate data generated using crystallographic or spectroscopy techniques. 
Computer programs used to generate such three dimensional models and/or 
perform the necessary fitting analyses include, but are not limited to: GRID (Oxford 
University, Oxford, UK), MCSS (Molecular Simulations, San Diego, CA), 
AUTODOCK (Scripps Research Institute, La Jolla, CA), DOCK (University of 
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California, San Francisco, CA), Flo99 (Thistlesoft, Morris Township, NJ), Ludi 
(Molecular Simulations, San Diego, CA), QUANTA (Molecular Simulations, San 
Diego, CA), Insight (Molecular Simulations, San Diego, CA), SYBYL (TRIPOS, Inc., 
St. Louis. MO) and LEAPFROG (TRIPOS, Inc., St. Louis, MO). 

5 In a preferred method of the present invention, the active site 

preferably comprises amino acid residues K343, E344, R347, R348 and D353 (or 
conservative substitutions thereof) according to Figure 6, amino acid residues E369 
and Y373 (or conservative substitutions thereof) according to Figure 6, amino acid 
residues E335, E386, E390, D398, E406, D407, E409 and E410 (or conservative 

10 substitutions thereof) according to Figure 6, or amino acid residues R358, R365, 
R368, R379, R380, R381, R384, R385, R394 and R397 (or conservative 
substitutions thereof) according to Figure 6. 

In the more preferred embodiment, the method of the present invention 
includes the use of an active site characterized by the three dimensional structure 

15 comprising the relative structural coordinates of amino acid residues K343, E344, 
R347, R348 and D353 according to Figure 8, ± a root mean square deviation from 

o 

the conserved backbone atoms of said amino acids of not more than 1.5 A, 

o 9 

preferably not more than 1.0 A, and most preferably not more than 0.5 A. In 
another embodiment, the active site is characterized by the three dimensional 
20 structure comprising the relative structural coordinates of amino acid residues 
E369 and Y373 according to Figure 8, ± a root mean square deviation from the 

o 

conserved backbone atoms of said amino acids of not more than 1.5 A, preferably 
not more than 1.0 A, and most preferably not more than 0.5 A. In another 
embodiment, the active site is characterized by the three dimensional structure 
25 comprising the relative structural coordinates of amino acid residues E335, E386, 
E390, D398, E406, D407, E409 and E410 according to Figure 8, ± a root mean 
square deviation from the conserved backbone atoms of said amino acids of not 
more than 1.5 A, preferably not more than 1.0 A, and most preferably not more 
than 0.5 A. In yet another embodiment, the active site is characterized by the 
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three dimensional structure comprising the relative structural coordinates of amino 
acid residues R358, R365, R368, R379, R380, R381, R384, R385, R394 and R397 
according to Figure 8, ± a root mean square deviation from the conserved 
backbone atoms of said amino acids of not more than 1.5 A, preferably not more 

5 than 1.0 A, and most preferably not more than 0.5 A. It is understood that the 
method of the present invention includes additional embodiments comprising 
conservative substitutions of the noted amino acids which result in the same 
structural coordinates of the corresponding residues in Figure 8 within the stated 
root mean square deviation. 

10 The effect of such an agent identified by computer fitting analyses on 

TNFR-1 DD, TNFR-1 DD complex or TNFR-1 DD analogue activity may be further 
evaluated computationally, or experimentally by competitive binding experiments 
or by contacting the identified agent with TNFR-1 DD (or a TNFR-1 DD complex or 
analogue or a molecule comprising TNFR-1 DD such as TNFR-1) and measuring the 

15 effect of the agent on the target's biological activity. Standard enzymatic assays 
may be performed and the results analyzed to determine whether the agent is an 
inhibitor of TNFR-1 DD activity (i.e., the agent may reduce or prevent binding 
affinity between TNFR-1 DD and a relevant binding protein, such as TRADD DD. 
Further tests may be performed to evaluate the selectivity of the identified agent to 

20 TNFR-1 DD with regard to other TNFR- 1 DD analogues or TRADD DD binding 
targets. 

An agent designed or selected to interact with TNFR-1 DD preferably is 
capable of both physically and structurally associating with TNFR-1 DD via various 
covalent and/or non-covalent molecular interactions, and of assuming a three 
25 dimensional configuration and orientation that complements the relevant active 

site of TNFR-1 DD. 

Accordingly, using these criteria, the structural coordinates of the 
TNFR-1 DD molecule as disclosed herein, and/or structural coordinates derived 
therefrom using molecular replacement or homology modeling, agents may be 
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designed to increase either or both of the potency and selectivity of known 
inhibitors, either by modifying the structure of known inhibitors or by designing 
new agents de novo via computational inspection of the three dimensional 
configuration and electrostatic potential of a TNFR-1 DD active site. 
5 Accordingly, in one embodiment of the invention, the structural 

coordinates of Figure 8 of the present invention, or structural coordinates derived 
therefrom using molecular replacement or homology modeling techniques as 
discussed above, are used to screen a database for agents that may act as potential 
inhibitors of TNFR-1 DD activity. Specifically, the obtained structural coordinates 
10 of the present invention are read into a software package and the three 
dimensional structure is analyzed graphically. A number of computational 
software packages may be used for the analysis of structural coordinates, including, 
but not limited to, Sybyl (Tripos Associates), QUANTA and XPLOR (Brunger, A.T., 
(1994) ggw 3_S5V a gygtpm for X-rav Crystallography and NMR. Xplor Version 
15 3.851 New Haven, Connecticut: Yale University Press) . Additional software 

programs check for the correctness of the coordinates with regard to features such 
as bond and atom types. If necessary, the three dimensional structure is modified 
and then energy minimized using the appropriate software until all of the structural 
parameters are at their equilibrium/optimal values. The energy minimized 
20 structure is superimposed against the original structure to make sure there are no 
significant deviations between the original and the energy minimized coordinates. 

The energy minimized coordinates of TNFR-1 DD bound to a "solved" 
inhibitor are then analyzed and the interactions between the solved ligand and 
TNFR-1 DD are identified. The final TNFR-1 DD structure is modified by 
25 graphically removing the solved inhibitor so that only TNFR-1 DD and a few 

residues of the solved agent are left for analysis of the binding site cavity. QSAR 
and SAR analysis and/or conformational analysis may be carried out to determine 
how other inhibitors compare to the solved inhibitor. The solved agent may be 
docked into the uncomplexed structure's binding site to be used as a template for 
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data base searching, using software to create excluded volume and distance 
restrained queries for the searches. Structures qualifying as hits are then screened 
for activity using standard assays and other methods known in the art. 

Further, once the specific interaction is determined between the solved 
inhibitor, docking studies with different inhibitors allow for the generation of initial 
models of new inhibitors bound to TNFR-1 DD. The integrity of these new models 
may be evaluated a number of ways, including constrained conformational analysis 
using molecular dynamics methods (i.e., where both TNFR-1 DD and the bound 
inhibitor are allowed to sample different three dimensional conformational states 
until the most favorable state is reached or found to exist between the protein and 
the bound agent). The final structure as proposed by the molecular dynamics 
analysis is analyzed visually to make sure that the model is in accord with known 
experimental SAR based on measured binding affinities. Once models are obtained 
of the original solved agent bound to TNFR-1 DD and computer models of other 
molecules bound to TNFR-1 DD, strategies are determined for designing 
modifications into the inhibitors to improve their activity and/or enhance their 
selectivity. 

Once a TNFR-1 DD binding agent has been optimally selected or 
designed, as described above, substitutions may then be made in some of its atoms 
or side groups in order to improve or modify its selectivity and binding properties. 
Generally, initial substitutions are conservative, i.e., the replacement group will 
have approximately the same size, shape, hydrophobicity and charge as the original 
group. Such substituted chemical compounds may then be analyzed for efficiency 
of fit TNFR-1 DD by the same computer methods described in detail above. 

Various molecular analysis and rational drug design techniques are 
further disclosed in U.S. Patent Nos. 5,834,228, 5,939,528 and 5,865,116, as well 
as in PCT Application No. PCT/US98/16879, published as WO 99/09148, the 
contents of which are hereby incorporated by reference. 
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The present invention may be better understood by reference to the 
following non-limiting Examples. The following Examples are presented in order to 
more fully illustrate the preferred embodiments of the invention, and should in no 
way be construed as limiting the scope of the present invention. 

5 

Example 1 

Tumor Necrosis Factor Receptor 1 (TNFR-1) death domain (DD) is the 
intracellular functional domain responsible for the receptor signaling activities. To 
understand the transduction mechanism of TNFR-1 signaling, we have performed 

10 structural and functional analysis of the TNFR-DD. The secondary structure of the 
TNFR-DD shows that it consists of six antiparallel helices. The determination of the 
topological fold and an extensive mutagenesis analysis revealed that there are two 
opposite faces that are involved in self-association and interaction with the TRADD 
death domain. Interestingly, the same critical residues in TNFR-DD are involved in 

15 both these interactions. There is a good correlation between the binding activities 
of the mutant proteins and their cytotoxic activities. These results provide 
important insight into the molecular interactions mediating TNFR-DD self- 
association and subsequent recruitment of TRADD in the signaling activity of 

TNFR-1. 

20 

Experimental Methods 

TNFR-DD cloning and site-directed mutagenesis. The DNA sequence coding from 
amino acid A316 to R426 from hTNFR-1 was cloned in pRSETB (Invitrogen) after 
25 amplification by PCR. The primer at the 5'end introduces an initiation site with an 
JVdel site upstream the methionine, and the primer at the 3'end introduces a His tag 
after R426. The sequence was confirmed by sequencing. Single point mutations 
were introduced using the chameleon double-stranded site directed mutagenesis 
(Stratagene) and verified by sequencing of the coding region. 
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Preparation of R347A mutant for NMR studies. 20 grams of wet cell paste expressing 
TNFR-DD R347A were resuspended in 300 ml lysis buffer containing 50 mM 
TrisCl, pH 8.8, 50 mM NaCl, and protease inhibitor tablets (Boehringer Mannheim 
GmbH, Mannheim, Germany). Cells were lysed by two passages through a 
Microfluidizer (Microfluidics Corporation, Newton, MA). Cellular debris were 
removed by centrifugation at 15000 g for 30 min. The soluble extract was then 
loaded onto a 2x20 cm Ni-NTA Superflow (Qiagen) column pre-equilibrated with 
lysis buffer. The column was washed with large amount of lysis buffer until no 
protein was present in the flow-through as judged by Bradford protein assay. 
TNFR-DD R347A was eluted with a step elution of 400 mM Imidazole in lysis 
buffer. DTT was added immediately to the eluate to a final concentration of 10 
mM. Death domain containing fractions were pooled and dialysed against buffer A 
(50 mM Tris.Cl, pH 8.8, 50 Imidazole, 10 mM DTT) overnight, and applied onto a 
1 x 10 cm Toyopearl QAE-550C anion exchange column (TosoHaas, 
Montgomeryville, PA) pre-equilibrated with buffer A. Greater than 95% pure 
R347A mutant was eluted with a linear NaCl gradient. The fractions were 
subjected to SDS-PAGE analysis and good fractions were combined and 
concentrated before loading onto a Superdex 75 gel filtration column (Amersham 
Pharmacia Biotech, Piscataway, NJ) equilibrated with buffer containing 50 mM 
Tris.Cl, pH 8.8, 100 mM NaCl, 50 mM Imidazole, 10 mM DTT. R347A containing 
fractions were eluted, dialysed and concentrated for NMR studies. 

NMR spectroscopy. All NMR spectra were acquired at 25°C on Varian Unity+ 600 
spectrometer. A series of HSQC [24] related three dimensional triple-resonance 
experiments were recorded using uniformly 15 N- and 15 N/ 13 C-labeled protein in 
water with 5% D 2 0, including HNCO, (HB) CBCACONNH, HNCACB, 
C(CC)_TOCSY_(CO)NNH, H(CC)_TOCSY_(CO)NNH, and HAHB(CBCA) CONNH, 
and HN (CACB) HAHB [23,24] to make backbone, side-chain and sequential 
assignments. 3D-HNHA [38] and 2D HMQC-J [50] were carried out for 3 J NHa 
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coupling constant measurement. R347A is a proline rich protein that contains a 
total of nine proline residues. The HBCBCA(CO)_N_CAHA and HACAJM [28], 
which correlates sequential residues via their 15 N nuclei were used to assign proline 
residues and additional amide residues which were missing due to fast exchange 
5 with water at pH 8.8. Validation and additional residue assignments were made 
using both 15 N-edited NOESY and 13 C-edited NOESY [29,27] experiments. Data 
sets were processed and displayed on SGI work station using the programs 
NMRDraw and NMRPipe [30]. The programs PIPP and STAPP [32] were used for 
peak picking, data analysis and part of auto-sequential assignment. 

10 

TNFR-DD three-dimensional structural model The TNFR-DD sequence alignment 
with Fas-DD was determined using QUANTA98 (Molecular Simulations, Inc., San 
Diego, CA). An initial homology model for TNFR-DD based on its sequence 
alignment with Fas-DD was accomplished with MODELLER [54, 55]. The TNFR- 

15 DD homology model was further defined to maintain consistency with the observed 
secondary structure elements observed by NMR. Basically, the structure was 
subjected to a simulated annealing protocol (40) using the program XPLOR (42) 
where only regions of the protein that differed from the secondary structure 
observed by NMR were free to move. (|> and i|j dihedral constraints based on 

20 accepted values for standard secondary structure elements were applied to residues 
identified to be a -helical by NMR, otherwise no constraints were applied and the 
residue was free to move. The Fas-DD sequence and solution structure (PDB 
accession code 1DDF) were used as reported by Huang, et al. [12]. 

25 Protein preparation for binding assay. Wild-type and mutant recombinant TNFR-DD 
proteins were prepared from Escherichia Coli strain BL21 (DE3) pLysS (Stratagene) . 
When cell growth reached an OD 595 of -0.5 at room temperature, protein 
expression was induced by adding 0.3-mM isopropyl-l-thio-p-D-galactopyranoside 
for 4 h. Cells were lysed with the french press (1000 psi) in 50 mM Tris pH 7.5, 50 
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mM NaCl, 10 jig/ml leupeptin, 10 jig/ml aprotinin, 2 mM phenylmethylsulfonyl 
fluoride at 4°C. The lysate was centrifuged at 21000 g at 4°C for 30 minutes. The 
soluble mutants were then incubated for 1 hour with Ni-NTA resin (Qiagen) at 4°C. 
The resin was washed twice in 50 mM Tris pH 7.5, 50 mM NaCl and twice in 50 

5 mM Tris pH 7.5, 50 mM NaCl, 50 mM imidazole. The bound protein was eluted 
from the resin in 50 mM Tris pH 7.5, 50 mM NaCl, 250 mM imidazole and 10 mM 
DTT was added to the eluate. The wild-type protein and most of the mutants were 
in the pellet fraction after lysis. The pellet was then washed twice in 50 mM Tris 
pH 8.5, 50 mM NaCl and the death domain proteins were extracted in 50 mM 

10 NaOAc pH 4.0, 50 mM MgS04. 10 mM DTT was then added to the extract. 

Effect of ionic strength on TNFR-DD self-association and association with TRADD-DD. 
The experiments were carried out according to the protocol described in the next 
paragraph (TNFR-DD binding assay). Only, the binding buffer was supplemented 
15 with NaCl 5M to obtain the desired final concentrations of 140 mM, 300 mM, 600 
mM and 1.2 M. 

TNFR-DD binding assay. TNFR-DD binding to itself or to TRADD was determined 
by an enzyme-linked immunosorbent assay (ELISA) . Ninety six well plates were 

20 coated overnight at 4°C with 200 ng of MBP-TNFR-DD, MBP-TRADD or MBP- 
MADD-LZ (N-terminal region of MADD containing a leucine zipper motif) in 
coating buffer (25 mM MOPS pH 7.5, 150 mM NaCl, 10 mM DTT, 0.02% NaN 3 ). 
The plates were then incubated in blocking buffer (10 mM MOPS 7.5, 150 mM 
NaCl, 0.05% Tween 20, 0.1% Gelatin, 0.02% NaN 3 ) for one hour followed by a one 

25 hour incubation with the different TNFR-DD proteins in binding buffer (100 mM 
Tris 7.5, 140 mM NaCl, 0.1 mM EDTA, 0.2% Triton X-100, 10 mM DTT). The 
plates were washed 4 times in 10 mM KP0 4 pH 7.4, 0.05% Tween 20. The bound 
His tagged proteins were then incubated for 30 minutes in 10 mM MOPS 7.5, 
150mM NaCl, 0.05% Tween 20, 0.02% NaN 3 , anti-6x His mAb (1:1000), anti- 
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mouse Biotin (1:1000), Streptavidin-Alkaline Phosphatase (1:1000). The plates 
were washed 4 times in 10 raM KP0 4 pH 7.4, 0.05% Tween 20. Tropix CDP-Star 
/Sapphire II solution was used for the luminescent reaction and detection with the 
lumicount plate reader (Packard) . 

Results 

Structure determination. Structure determination of the wild-type death domain of 
TNFR-1 by NMR has been hindered due to low protein solubility over a wide range 
of pH (4-10). Therefore, a series of single point mutants have been made and 
tested for their solubility at different pH. The mutant R347A was found to be 
soluble and stable in concentrations up to 14 mg/ml (~1 mM) at pH 8.8. Under 
these conditions, the protein was shown to be predominandy monomeric by size 
exclusion chromatography. The HSQC spectrum of 15 N-labeled protein indicates 
that the protein was folded. Therefore, this mutant protein was chosen for the 

subsequent NMR studies. 

However, many of the amide proton peaks were severely broadened 
due to rapid exchange with water at high pH (pH 8.8) which hindered the 
complete assignment of TNFR-DD by HSQC-related triple resonance experiments. 
Nevertheless, the majority of TNFR-DD forms a-helices, where the amide proton 
exchange is reduced by hydrogen bond formation, allowing for the observation of 
backbone chemical shifts and the determination of sequential assignments in the 
helical regions. A total of 63 residues were assigned in the region from A328 to 
C413, where the amide signals of terminal residues, A316-P327 and G414-R426 
were not be observed. In an effort to complete the sequential assignment, two 
experiments were designed based on recently developed methodologies where 
NMR experiments correlated sequential residues via 15 N and Ha chemical shifts 
instead of the amide proton. These experiments enabled us to assign 24 additional 
residues. Therefore, a total of 87 residues out of 112 amino acids were assigned. 
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Figure 1 summarizes the NMR secondary structure data for TNFR-DD 
R347A. This data indicates that TNFR-DD consists of six a-helices where helices 1, 
2, 4, 5 and 6 exhibit the typical NOE pattern for an helix ( 3 JHN K < 5Hz, short- and 
medium-range NOEs: d NN (i,i+l), d aN (i,i+3) d aN (i,i+4) and d ap (i,i+3) and Ca/CP 

5 secondary chemical shifts) . These helices are separated by loop regions that were 
not observed because of the lack of amide proton resonance peaks due to the rapid 
solvent exchange. Ca/Cp chemical shift deviation indicates that both terminal 
regions are unstructured. A comparison of the secondary structure of TNFR-DD 
with Fas-DD indicates that the main structural difference between the two proteins 

10 resides in the length of helix-3. Helix 3 in the Fas-DD structure encompass nine 
residues from D240 to E248, whereas in TNFR-DD R347A, only three amide 
protons (E355 to D357) were detectable. The lack of amide peaks in the C- 
terminus region hinders the assignment of the secondary structure for these 
residues. However, the relatively high amide exchange rate suggests that the 

15 structure of this region is somewhat disordered, possibly due to a higher flexibility 
of helix-3. Data obtained by D 2 0 exchange with the mutant R347K at pH4 supports 
this observation (data not shown). Interestingly, similar results were obtained with 
FADD-DD [13]. The two dimensional helical wheel topology shown in Figure 2 is 

1 3 

based on more than 40 long range inter-helical side chain NOEs assigned from C 
20 edited NOESY experiments, combined with secondary structure analysis (Fig. 1) . 
TNFR-DD is consistent with other published DD structures where the protein 
consists of six helices arranged in an anti-parallel fashion. Many of the long range 
NOEs were found between side chain residues of a2 making contacts with side 
chain residues of «3 (K343-I356, V346-I356, A347-I356), CC4 (V346-M374, L349- 
25 M374, L349-W378) and CC5 (F345-L389, L349-L392, L349-V395, R348-V395). A 
few NOE pairs involving contacts between al and a4 (V334-Q371) as well as 
between al and a6 (V333-I408 and V337-I408) were also found. Only one side 
chain contact between a5 and a6 was observed, with L401 at the C-terminus of «5 
and L405 at the N-terminus of a6, suggesting that the angle between these two 
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helices is such that it prevents additional side chain contacts between them. Based 
on the planar view centered on helix-2 (Fig. 2), helices 2-5 form a twisted four 
helix bundle with a hydrophobic core consisting of residues F345, V346, 1356, 
M374, L392. The side chain of the hydrophilic residues in the helical bundle are 
solvent exposed. The long range inter-helical NOE between side chain of the 
residues in the core region indicates that these residues are close to each other and 
probably form hydrophobic interactions to stabilize the protein folding. 

The secondary structure assignment of R347A is very similar to that of 
the mutant R347K, which has been shown to have a similar structure to that of the 
wild-type protein, based on HSQC spectra comparison performed at pH4.0 (data 
not shown) . Therefore, the overall structure of R347A is most likely similar to that 
of the wild-type. 

Effect of ionic strength on TNFR-DD self-association and its interaction with TRADD- 
DD. Based on the observation that charged residues have been shown to be 
involved in the Fas-FADD death domains interaction, we examined the effect of 
sodium chloride on the TNFR-DD self association and the interaction with TRADD 
death domains. Figure 3 shows that increasing concentrations of NaCl are 
accompanied by a decrease in death domain interaction in both cases, implying 
that electrostatic interactions and/or hydrogen bonds are important for the TNFR- 
death domain self-association and its interaction with TRADD-DD. 

Effect of TNFR-DD single point mutations on self-association. Since electrostatic 
interactions might play an important role in death domain interactions, a series of 
19 charged amino acids were mutagenized (Fig. 4). Alanine was chosen as the 
substituent amino acid, since it eliminates the contribution from the side chain of 
the replaced amino acid with a potential minimal alteration of the overall protein 
structure. For tyrosine 372, the substitution was made with a phenylalanine in 
order to abolish the contribution of the hydroxyl group while keeping the aromatic 
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side chain intact. The wild-type and mutant constructs were expressed in bacteria 
as His-tagged proteins at the C-terminus. The wild-type and most of the mutant 
proteins were in the insoluble fraction after lysis at pH7.5. Only mutants E344A 
and R347A, were soluble upon lysis at pH7.5. These mutant proteins were purified 

5 from the lysate with Ni-NTA resin. All the others including the wild-type were 
extracted directly from the pellet fraction at pH 4.0. The purity of the proteins was 
determined to be about 80% by SDS-PAGE. 

All the mutants depicted in Figure 4 were studied for their ability to 
interact with the wild-type death domain in a cell free binding assay. The wild-type 

10 TNFR-DD was expressed as a MBP fusion protein (MBP-TNFR-DD) and was 

immobilized on 96 well plates for ELISA. Despite the tendency of the wild type and 
most of the mutant proteins to aggregate at pH 7.5, all of them were soluble at that 
pH over the concentration range that was used for these binding experiments (0.1- 
5 |lg/ml) . The four mutations K343A, E344A, R347A and R348A in CC2 abolish the 

15 interaction with the wild-type death domain (Fig. 5A). Mutation of the first residue 
of Ct3, D353A, also abolishes the interaction with the wild-type death domain. 
Conversely, the next three sequential mutations in <X3, D357A, R358A and E360A, 
and R365A in loop 3 retain an interaction with MBP-TNFR-DD similar to that of the 
wild-type protein. The first two sequential mutations in cc4, E369A and Y372F, lost 

20 their interaction with the wild-type death domain similarly to the mutants in a2. 
All the other mutations in 0C4 (R379A, R380A, R381A), Loop4 (R384A, E386A), (X5 
(E390A) and a6 (E406A, E410A) retained an interaction with MBP-TNFR-DD 
similar to that of ±e wild-type death domain. 

25 Effect of TNFR-DD single point mutations on interaction with TRADD-DD. The TNFR- 
DD mutants were also studied for their ability to interact with TRADD. The full 
length TRADD protein was expressed as a MBP fusion protein (MBP-TRADD) and 
was immobilized on 96 well plates for ELISA. The trend of the binding ability of all 
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the different mutants to MBP-TRADD is the same as for the binding to MBP-TNFR- 
DD (Fig. 5B). 



Discussion 

5 We have performed NMR studies and extensive mutagenesis analysis of 

the TNFR-DD. Overall, the secondary structure of TNFR-DD is similar to the 
structure of other death domains. It consists of six antiparallel a-helices. 
Interestingly, it appears that a2, part of a3 and part of 0t4 contain important 
residues for the self-association and interaction with TRADD-DD. Individual 
10 mutation of the four charged residues in a2, K343A, E344A, R347A and R348A 
results in a loss of interaction with TNFR-DD and with TRADD-DD. These residues 
are solvent exposed and located on the same side of the protein according to the 
deduced topological fold (Fig. 2). On the other hand, only the first mutated 
residue in a3, D353A, has the same effect. The other three mutations in 063 
15 (D357A, R358A and E360A) have no or little effect on the binding properties of the 
TNFR-DD. In the case of a4, the first two mutants E369A and Y372F are deficient 
in their interaction with TNFR- death domain and TRADD-DD. The other three 
mutations, R379A, R380A, R381A retain binding properties that are similar to the 
wild-type protein. The other mutants in loop 3 (R365A), loop 4 (R384A and 
20 E386A), a5 (E390A) and Ct6 (E406A and E410A) also seem to be unaffected in 
their binding properties when compared to the wild-type protein. These results are 
in good agreement with the cytotoxicity data obtained previously by [7, 10]. 
Mutants R347A and E369A that are defective in our ELISA binding assays are also 
deficient in cytotoxicity. In contrast, mutants R379A, E390A, E406A and E410A 
25 that have unaltered binding properties also retained cytotoxic activity. Only 

mutant K343A that lost its ability to interact with the wild-type receptor DD or with 
TRADD-DD retained some cytotoxic activity, also to a lesser extent than the wild- 
type receptor [7, 10] . The correlation between death domain interaction and 
cytotoxicity further supports the hypothesis that the aggregation of TNFR- 1 and 
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subsequent recruitment of TRADD through death domain interactions is critical for 
generating a cytotoxic signal [8, 51, 52]. From our results, it is possible that self- 
association and binding to TRADD-DD use similar surfaces since all the mutations 
have the same effect in both cases. When these mutations are placed in the context 
5 of the topological fold, it appears that there are two opposite and solvent exposed 
faces involved in the interaction. On one side there are critical residues in a2 
(K343A, E344A, R347A and R348A) and the beginning of (X3 (D353A) and on the 
other side there are critical residues in ct4 (E369A, Y372F). It was recently shown 
that a2 and a3 of FADD-DD contain critical residues for the interaction with Fas- 
10 DD [13]. However, the role of a4 in FADD-DD has yet to be addressed. 
3 Interestingly, a comparison of mutations in a3 of Fas and TNFR-1 death 

0 domains reveals differences in the contribution of this helix to their respective 

p binding properties. Regarding the self-association, the four mutations in 0C3 of Fas 

| (E240A, D244A, E245A and K247A) abrogate their self-association [12] . On the 

n 15 contrary, in the case of TNFR-1, for the four equivalent mutations, only the first 
3 one in oc3 (D353A) abrogates the interaction with the wild type death domain, as 

2 the next three mutations (D357A, R358A and E360A) do not (Fig. 4). Therefore, it 

t would seem that the entire helix-3 of Fas-DD is important for the self-association, 

* ' as in the case of TNFR-DD only the beginning of helix-3 is important for the self- 

20 association. In the case of the Fas-FADD interaction, the first two mutations in 0,3 
of Fas-DD, E240A and D244A gready affect their interaction with FADD-DD, as the 
next two mutations E245A and K247A only have a marginal effect. In the case of 
TNFR-1-TRADD interaction, only the first mutation in a3 of TNFR-DD (D353A) 
abrogates the interaction with TRADD-DD, as the next three mutations (D357A, 
25 R358A and E360A) do not. It is worth mentioning that while self-association of 
Fas-DD may not be required for the recruitment of FADD to the receptor, as 
suggested by the mutants E245A and K247A, we were unable to find mutations in 
TNFR-DD that have a differential effect on self-association and TRADD binding. 
Thus, it appears that a3 contains critical residues that can account for the 
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differences between TNFR-1 and Fas regarding their self-association and 
interaction with other death domain proteins. 

The difference between death domain interaction pairs is further 
demonstrated by the complex of the death domains of Pelle and Tube. The three 

5 dimensional structure of a complex between the death domains of Pelle and Tube 
was recently published [19]. It shows that the surfaces involved in the interaction 
seem to be different between different pairs of death domains. The Tube-Pelle 
dimers rely on contacts between a4 and a5 of Pelle with a6 and a C-terminal tail 
of Tube. This is in contrast with the homo- and hetero-association of the death 

10 domains of TNFR-1 and TRADD, which mainly rely on a2, 063 and a4 of TNFR-1 . 
These results suggest that while a topological fold is maintained between members 
of the death domain family, the observed selectivity and specificity is obtain by a 
combination of distinct sequences and surface interpositions between the death 
domains. 

15 Because of the charged nature of the death domains and the caspase 

recruitment domains, it was previously proposed that electrostatic interactions 
might mediate the homotypic interaction [12, 13]. However, the crystal structure 
of the CARD of Apaf-1 in complex with the CARD of procaspase-9 has been solved, 
showing that the specificity as well as the driving force for this interaction is 

20 provided by a combination of extensive hydrogen bond contacts and van der Waals 
interactions [18]. Moreover, this interaction is insensitive to high ionic strength, 
confirming that the driving force of the interaction is not mediated by electrostatic 
interactions. Similarly in the case of Pelle and Tube death domains interaction, it 
has been shown that hydrogen bonds and van der Waals interactions are the 

25 driving force for this interaction [19] . In the case of TNFR-DD, we have shown that 
interaction with the wild-type death domain and interaction with TRADD-DD is 
disrupted at high salt concentration (0.3 M-1.2M). These results suggest that 
TNFR-DD interactions are mediated, at least in part, by electrostatic interactions, 
similarly to what has been shown in the case of Fas-DD [12]. This is also 
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supported by our mutagenesis study in which mutation of many of the charged 
residues to alanine induces a loss of self interaction or interaction with TRADD-DD. 
These different studies show that not only the surfaces, but also the nature of the 
interactions involved in death domain associations is not conserved, which provides 

5 for a mechanism of death domain selectivity and specificity of interaction. 

A homology model for TNFR-DD was built utilizing the Fas-DD three- 
dimensional structure [12] . This model of TNFR-DD as well as the two 
dimensional helical wheel topology (Fig. 2) show that there are two opposite faces 
involved in TNFR-DD self-association and interaction with TRADD-DD. In TNFR-1 

10 signaling, the receptor trimerizes upon binding with a TNF trimer and then recruits 
TRADD. Our experiments show that the critical residues for TNFR-DD self- 
association or interaction with TRADD-DD are the same. Therefore, it is possible 
that self-association of TNFR-DD and interaction with TRADD-DD occurs via the 
same binding surface. If this is the case, it is possible that ligand-induced 

15 trimerization of the receptor recruits TRADD by simply increasing the affinity for 
TRADD due to an effective increase in the local concentration of TNFR-DD. 
Alternatively, TNFR-1 trimerization may create a high affinity-binding site for 
TRADD. It is also possible that a third partner is involved in the regulation of the 
recruitment of TRADD to TNFR-1; Recently, SODD was shown to be constitutively 

20 bound to TNFR-1 intracellular domain and released from the receptor upon binding 
of the ligand, allowing for TRADD to be recruited to TNFR-1 [53]. It is possible 
that SODD prevents TRADD from interacting with TNFR-1 in the absence of ligand 

bound to the receptor. 

In conclusion, TNFR-DD consists of six antiparallel helices similar to 
25 other death domains where the major structural and functional difference occurs in 
CC3. Based on mutagenesis study, helices a2 and part of oc3 constitute a binding 
surface on one side of the protein, while CC4 constitutes another binding surface on 
the opposite side of the protein. The critical residues involved in self-association as 
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well as interaction with TRADD-DD are the same, suggesting a very similar mode of 
interaction in both cases. 

Example 2 

5 The solution structure of the R347K mutant of the TNFR-1 DD was 

solved by NMR spectroscopy. A total of 20 structures were calculated by means of 
hybrid distance geometry-simulated annealing using a total of 1360 experimental 
NMR restraints. The atomic rms distribution about the mean coordinate positions 
for the 20 structures for residues composing the six a-helices is 0.40 A for the 

10 backbone atoms and 1.09 A for all atoms. The structure consists of six anti-parallel 
a-helices arranged in a similar fashion to the other members of the death domain 
superfamily. The secondary structure of R347K TNFR1-DD is very similar to the 
secondary structure deduced for the R347A TNFR1-DD mutant (Example 1) with 
the exception of the length of helix 3. Helix 3 was observed to be shorter by three 

15 residues in the R347A TNFR1-DD mutant presumably a result of the inherent 
flexibility of this helix and the increase amide exchange rate at pH 8.8. The 3- 
dimensional structure of R347K TNFR-1 DD is consistent with the deduced 
topology of R347A TNFR-1 DD. Mutagenesis studies set forth in Example 1 
identified critical residues located in helix 2 and part of helices 3 and 4 that are 

20 crucial for self-interaction and interaction with the TRADD DD. Structural 
superposition with previously solved proteins in the death domain superfamily 
reveals that the major differences between the structures reside in helices 2, 3, and 
4. 

25 Experimental Methods 

TNFR-DD cloning and site-directed mutagenesis (R347K TNFR-1 DD). The DNA 
sequence coding from amino acid 316 to 426 from hTNFR-1 was cloned in pRSETB 
(Invitrogen) after amplification by PCR. The primer at the 5'end introduces an 
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initiation site with an Ndel site upstream of the methionine, and the primer at the 
3'end introduces a His tag after amino acid 426. The sequence was confirmed by 
sequencing analysis. A single point mutation (R347K) were introduced using the 
chameleon double-stranded site directed mutagenesis (Stratagene) and verified by 
5 sequencing of the coding region. 

Protein expression. R347K TNFR-1 DD protein was overexpressed in BL-21 E. coli. 
Cells were grown at 25°C on minimal medium containing either 2g/L 15 N 
ammonium sulfate or 15 N ammonium sulfate and 2 g/L 13 C glucose. The cells were 
10 then lysed and R347K TNFR-1 DD mutant protein was purified (see purification 
section) . 

Purification ofR347K TNFR-1 DD protein. 20 grams of wet cell paste expressing 
R347K TNFR-1 DD were resuspended in 300 ml lysis buffer containing 50 mM 

15 TrisCl, pH 7.5, 50 mM NaCl, and protease inhibitor tablets (Boehringer Mnnheim 
GmbH, Mannheim, Germany). Cells were lysed by two passages through a 
Microfluidizer (Microfluidics Corporation, Newton, MA). The cell lysate was 
centrifuged at 15000 g for 30 min. The TNFR-1 DD protein was then extracted 
from the pellet with 200 mL pH 8.8 buffer (50mM Tris pH 8.8 and 50 mM NaCl). 

20 The extraction was performed three times until no additional protein was 

solubilized. Forty mL of Ni-NTA resin (Qiagen) was then added to the combined 
extract. The resin was first washed with pH 8.8 buffer followed by pH 8.8 buffer 
containing 50 mM imidazole. The R347K TNFR DD protein was then eluted with 
400 mM imidazole in pH 8.8 buffer. DTT was added immediately to the eluate to a 

25 final concentration of 10 mM. The eluted protein was pooled and dialyzed against 
buffer containing 50 mM sodium acetate pH 4.0-4.2 and 10 mM DTT overnight. 
The protein was then concentrated using centriprep (YM-10, Amicon) to about 6 
mg/ml and applied to G3000SW (TOSOHASS) column pre-equilibrated with lOmM 
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sodium acetate pH 4.0 and 10 mM DTT. The monomeric fraction was pooled and 
concentrated for NMR studies. 



NMR data collection, complete assignments and secondary structure determination. All 

5 spectra were recorded at 35° C on a Varian Unity Plus 600 spectrometer equipped 
with a triple-resonance ^/^C/^N probe and an actively shielded z-gradient pulsed 
field accessory. X H- 15 N HSQC and all HSQC-based 15 N edited 3D experiments on 15 N 
labeled and triple-resonance 3D experiments on 15 N/ 13 C double-labeled protein 
were recorded with the enhanced-sensitivity pulsed field gradient approach [23] . 

10 This approach provides coherence transfer selection both to improve sensitivity and 
to eliminate artifacts as well as for solvent suppression. The complete assignments 
(>95%) of the 1 H, 15 N and 13 C resonances were based on the following 
experiments: CBCA(CO)NNH, HNCACB, C(CC)TOCSY_NNH, H(CC)TOCSY_NNH, 
HAHB (CBC ACO) NNH [23, 24] [Table IS]. 2D 13 C (methyl) HSQC and methyl 

15 relay experiments were used for auxiliary methyl assignments of valine, leucine and 
isoleucine residues [25, 26]. Some ambiguous resonances were further confirmed 
by analysis of the simultaneous 15 N/ 13 C-edited NOESY [27] experiment. R347K 
TNFR-1 DD is a proline rich protein, which contains several peptidyl-prolyl 
(Cis/Trans) isomers. Most of the prolines were detected using a new experiment 

20 that connects sequential residues via their 15 N nuclei [28] . 4D 13 C, 13 C -edited 
NOESY [29] were carried out on a double labeled ( 15 N/ 13 C) sample in D 2 0 to 
assign overlapped y -methyls of valine and the 6 -methyls of leucine and isoleucine 
residues and also to identify some long-range NOEs between these methyls. Data 
sets were typically processed and displayed on an SGI workstation using the 

25 NMRDraw and NMRPipe [30] programs. A skewed 60° phase-shifted sine-bell 

function and a single zero-filling were used in each of the three dimensions prior to 
Fourier transformation. For triple-resonance 3D experiments, the time domain was 
extended by a factor of two using froward-backward linear prediction in the N 
(t2) dimension. For constant-time ^-^C correlation experiments, mirror image 
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linear prediction was used prior to zero-filling to double the time domain data 
points [31]. The programs PIPP and STAPP [32] were used for data analysis and 
semi-automatic assignments [31]. 

The secondary structure of R347K TNFR-1 DD was based on 

5 characteristic NOEs involving the amide NH, H„ and Hp protons from 15 N edited 
NOESY and simultaneous 15 N/ 13 C-edited NOESY spectra, three-bond J-coupling 
constants 3 J NHct from the 3D HNHA, slowly exchanging amide protons with D 2 0, 
and C K and C p secondary chemical shifts. It was determined that the secondary 
structure of R347K TNFR-1 DD is composed of six helices, ranging in length from 7- 

10 10 residues. Helix 3 (residues D353-L361) appears to be more flexible as evident 
by the increased amide exchange rates relative to other helical regions in the 
protein. 

R347K TNFR-1 DD Structure Calculation. The NMR solution structure is based on 
15 inter-proton distance constraints converted from observed NOEs in both the 15 N- 
edited NOESY [33, 34] and simultaneous 15 N/ 13 C-edited NOESY experiments. The 
NOEs were classified as either strong (1.8-2.7 A), medium (1.8-3.3 A) or weak 
(1.8-5.5 A) constraints. Upper distance limits for distances involving methyl 
protons and non-stereospecifically assigned methylene protons were corrected 
20 appropriately for center averaging [35], and an additional 0.5 A was added to 

upper distance limits for NOEs involving methyl protons [36, 37]. <j) and i|f torsion 
angle constraints were obtained from the 3D HNHA experiment [38] and from 15 N, 
Hot, C a and C p chemical shifts using the TALOS [39] program. 

The structures were calculated using the hybrid distance geometry- 
25 dynamical simulated annealing method of Nilges et al. [40] with minor 
modifications [41] using the program XPLOR [42], adapted to incorporate 
pseudopotential secondary C a /Cp chemical shift restraints [43] and a 
conformational database potential [44, 45] . The target function that is minimized 
during restrained minimization and simulated annealing comprises only quadratic 
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harmonic terms for covalent geometry and secondary C«/Cp chemical shift 
restraints, square-well quadratic potentials for the experimental distance and 
torsion angle restraints, and a quartic van der Waals term for non-bonded contacts. 
All peptide bonds were constrained to be planar and trans. There were no 
hydrogen-bonding, electrostatic or 6-12 Lennard-Jones empirical potential energy 
terms in the target function. 

The N- and C-terminal regions of the protein were unstructured and 
contain several peptidyl-prolyl (Cis/Trans) isomers therefore, they were not used in 
the structure calculations. The structure of R347K TNFR-1 DD encompassing 
residues P327-D413 was determined from a total of 1167 distance constraints 
comprising 264 intra-residue, 279 sequential, 247 medium, and 295 long range 
distance constraints, 82 hydrogen bond constraints and 112 torsion angles 
constrains comprised of 56 (|) and 56 ijj dihedral constraints. The hydrogen bond 
constraints were based on the observation of slow exchanging NH protons in a D 2 0 
solution monitored by an J H- 15 N HSQC spectrum. 

The final ensemble of 20 structures contained no distance constraint 
violations greater than 0.25 A and no torsion angle constraint violations greater 
than 3°. The NMR structures are well defined, as evident by the atomic r.m.s.d. of 
the 20 simulated annealing structures about the mean coordinate positions where 
the backbone of the helical regions and all atom in the helical region is 0.40 A and 
1.09 A, respectively (Table 1). The R347K TNFR-1 DD NMR structure is consistent 
with a good quality structure based on PROCHECK and Ramachandran analysis 
[46, 47] . A Ramachandran plot of the minimized average structure shows a total 
of 88.5% of the residues are in the most favored region and 11.5 % in the 
additional allowed region. PROCHECK analysis indicates an overall G-factor of 
0.01 and only 1 bad contacts. 
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Results and Discussion 

Structure Determination. Historically, death domains have posed problems for NMR 
studies at physiological pH due to their tendency to self-associate and form large 
molecular weight aggregates. This is also the case for the Death Effector Domains 
(DED) and the Caspase Recruiting Domains (CARD) which are structurally related 
to the death domains. Therefore, structural studies of the death domain 
superfamily of proteins requires relatively low (< 4) or high O 8) pH to minimize 
the natural tendency for the proteins to self-aggregate. In addition, single point 
mutants which disrupt self-association were necessary for the structural studies of 
the FADD DED [16] and in the current study of TNFR-1 DD. 

The wild type TNFR-1 DD was very insoluble and not amenable to NMR 
structural studies over a wide pH range (pH 4-10), therefore mutant proteins were 
made to improve the solubility and stability of the protein. As previously described, 
the solubility and stability of the R347A mutant at high pH (pH=8.8) allowed for 
secondary structure determination and an insight into the topological fold of the 
protein (see. Example 1) . Due to severe broadening of the amide proton peaks 
because of rapid exchange with water at high pH (pH=8.8), complete assignment 
and structure determination of the protein was hindered. Nevertheless, it was 
demonstrated that the majority of TNFR-DD forms (X-helices, where the amide 
proton exchange is reduced by hydrogen bond formation. This allowed for the 
observation of backbone chemical shifts and the determination of sequential 
assignments in the helical regions using triple resonance experiments. The 
mutation of R347 to lysine allowed for NMR studies to be performed at low pH 
(pH=4.0) where the exchange with water is reduced. At pH=4.0, the solubility of 
the R347K TNFR-1 DD mutant was ~1 mM and the protein ran predominantly as a 
monomer on a gel filtration column. Additionally, the protein appears folded by 
HSQC correlation experiment collected on a 15 N labeled protein. 
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Figure 6 is a summary of NMR derived data that summarize the 
secondary structural elements of R347K TNFR-1 DD. Sequential and medium and 
range NOE's: d NN (i,i+l), d KN (i,i+3) d aN (i,i+4) and d ap (i,i+3) along with aC„ and 
aCr secondary chemical shifts identify the 6 helical regions of the protein. 
5 Additionally, the experimentally measured coupling constants ( 3 J HN h«) and tne 
hydrogen-deuterium exchange data correlate well with the 6 a-helices with the 
exception of helix 3 which shows that the amides located in helix 3 exchange 
quickly with D 2 0. The approximate lengths of the helices are from residues: A328 
to N336, W342 to L349, D353 to L361, L367 to R380, L389 to D398, and G403 to 
10 L412. The N- and C- terminal regions are unstructured and both exhibit multiple 
isomers (cis-trans) due to the presence of six proline residues. 

The structure of the R347K TNFR-1 DD is well defined by 1360 
experimentally derived NMR restraints (Table 1) consisting of 1167 NOE distance 
constraints, 112 dihedral angle constraints and 82 hydrogen bond constraints. An 
15 overlay of the backbone atoms for the best 20 simulated annealing structures for 
the R347K TNFR-1 DD was generated (not shown). As evident from the best-fit 
superposition of the R347K TNFR-1 DD NMR structures, the backbone and 
sidechains for the core of the protein consisting of the six helical regions are well 
defined having an r.m.s.d. of 0.40 A and 1.09 A, respectively. A ribbon diagram of 
20 the C K trace of the average minimized structure of R347K TNFR-1 DD was also 
generated (not shown) . The structure consists of six anti-parallel amphipathic a- 
helices where the majority of the hydrophobic residues form an extensive network 
of interactions that constitute the core of the protein. Helices 1,5, and 6 are 
aligned anti-parallel to one another and are aligned almost perpendicular to helices 
25 2, 3, and 4. The alignment of helix 3 is slightly skewed from being anti-parallel to 
helices 2 and 4. The surface of the R347K TNFR-1 DD mutant is composed of 
several basic (R341, K343, K347, R348, R358, R365, R368, R379, R380, R381, 
R384, R385, R394, and R397) and several acidic (E335, E344, D353, E355, D357, 
E360, E369, E386, E390, D398, D400, E406, D407, E409 and E410) residues. One 
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face of the R347K TNFR-1 DD surface contains an acidic patch consisting of 
residues from helix 1 (E335), helix 5 (E386, E390, D398) and helix 6 (E406, D407, 
E409, and E410). The other face of the protein contains several exposed basic 
residues from helix 4 (R365, R368, R379, R380, and R381), the loop between helix 
5 4 and helix 5 (R384 and R385), and from helix 5 (R394 and R397). The structure 
also contains several surface exposed hydrophobic residues (Y331, L340, L349, 
L351, L361, L367, Y372, A376, L391, M399, L402, and L412). These electrostatic 
patches on the surface of TNFR-1 DD is suggestive of potential protein binding 
sites. 

10 

Structural comparison of TNFR-1 DD mutants. Due to the poor solubility of the wild 
type TNFR-1 death domain, NMR structural studies were pursued on two different 
mutant proteins at two distinct pH values. The R347A mutant was studied at high 
pH (pH 8.8) where the protein is soluble to 1 mM (see, Example 1); whereas the 

15 R347K TNFR-1 DD described herein is soluble to 1 mM at pH=4.0. The mutation 
in both proteins occurs in helix 2 and although the solubilities of the R347A and 
R347K TNFR-1 DD proteins are drastically different, the secondary structure in this 
region appears unaffected. Additionally, the pattern of inter-helical long range 
NOEs is conserved when comparing the two TNFR-1 DD mutant proteins further 

20 supporting the observation that the structures of the two mutant proteins are very 
similar. The only notable difference between the R347A and R347K TNFR-1 DD 
proteins resides in the length of helix 3. The helix spans 5 residues (E353-D357) in 
the R347A mutant compared to 8 residues (D353-L361) in the R347K TNFR-1 DD. 
The observed difference in the length of helix 3 for the two mutants may be 

25 attributed to an apparent increase in the mobility of this helix relative to the other 
helical regions as suggested by the exchange data (Figure 6) . At the relatively high 
pH used for the analysis of the R347A mutant structure, the amide exchange rate 
with solvent is rapid compared to the conditions used for the R347K TNFR-1 DD 
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mutant. Therefore, it is not surprising that the remainder of helix 3 was not 
observed in the R347A study (see, Example 1). 

Structural comparison ofTNFR-1 DD with other members of the death domain 
superfamily. The death domain superfamily of proteins includes the death domains, 
the death effector domains, and the CARD domains which are involved in signaling 
programmed cell death and homotypic interactions. Each known structure for the 
proteins in this superfamily contains a characteristic core of six anti-parallel a- 
helices [12-14, 16-19]. Typically there is -20% sequence identity between the 
various death domains [48] and less sequence homology when comparing the 
death domains to the death effector domains and CARD domains. All of the 
members in this class contain relatively short protein sequences (-100 amino 
acids) that encompass the 6 a-helices. Figure 7 presents a pseudo-structural and 
sequence alignment for members of the death domain superfamily of proteins. 

A comparison of the R347K TNFR1 DD structure with the other known 
structures in the death domain superfamily has revealed that the overall fold for 
TNFR-1 DD is similar to other death domains. The superposition of the 6 a-helical 
regions of R347K TNFR-1 DD with the proteins in the death domain superfamily 
was generated (not shown) . The observed deviations when comparing the R347K 
TNFR-1 DD structure to the other protein structures in the death domain 
superfamily reside in either helix length (Figure 7) or helix orientation (not 
shown). A detailed analysis of the death domain structures revealed that the 
structure of p75 neurotrophin was the most similar to TNFR-1 DD having an 
r.m.s.d of 1.80 A for the 68 C a atoms composed of the helical core. The similarity 
to TNFR-1 DD is followed by the FAS DD and FADD DD having r.m.s.d's of 2.17 A 
and 2.14 A for the 70 C a and 52 C a atoms in the helical core, respectively. 
Comparison of TNFR-1 DD with the CARD domains resulted in r.m.s.d's of 1.95 A 
and 2.17 A for the 49 C a and 53 C a atoms in the core helical region for the RAIDD 
CARD and Apaf-1 structures, respectively. Overlay of R347K TNFR-1 DD mutant 
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with the FADD DED resulted in a r.m.s.d. of 2.08 A for the 56 C ft atoms in the 
helical regions of the proteins. 

Although the structural superpostition of R347K TNFR-1 DD with the 
other proteins in the death domain superfamily yielded reasonable r.m.s.d's 

5 consistent with the similarity in the overall folds, the major structural difference 
between R347K TNFR-1 DD and the other death domain-like proteins resides in the 
relative orientations of helices 2, 3, and 4. In the R347K TNFR-1 DD structure, 
helices 2 and 4 are aligned parallel to one another and the orientation of helix 3 is 
almost perpendicular. P75 neurotrophin is the only member of the death domain 

10 superfamily that retains the orientation of helices 2, 3, and 4. The lengths of 

helices 2 and 3 are very similar among all the family members, however, helix 4 of 
FAS DD [12], p75 neurotrophin [14], and FADD DED [16] are much shorter 
(consisting of 11, 10, and 7 residues respectively) than the length of helix for 
R347K TNFR1-DD. Additionally, the length of helix 1 of R347K TNFR-1 DD 

15 consisting of 8 residues is much shorter than helix 1 of the FADD DD (14 residues) 
[13], the RAIDD CARD domain (13 residues) [17], and Apaf-1 (12 residues) [21]. 

Several of the proteins in the death domain superfamily (FAS DD [12], 
RAIDD CARD [17], and the FADD DD [13]) have charged surfaces and it is 
presumed that electrostatic interactions are the driving force for homotypic 

20 interactions. Recently, crystal structures of complexes formed between procaspase- 
9 with APAF-1 CARD domains [18] and the Pelle DD with the Tube DD [19] have 
shown that the surfaces involved in the homotypic interactions are charged and of 
opposite polarity. Nevertheless, the mode of interaction in both cases was shown to 
be driven by a network of hydrogen bond contacts as well as by van der Waals 

25 interactions instead of an electrostatic interaction as suggested by the protein 
surface. In the case of the TNFR-1 DD, the interaction between TNFR-1 DD and 
TRADD DD could be disrupted by salt (see, Example 1). This is consistent with the 
molecular surface for the NMR structure of the R347K TNFR-1 DD mutant which 
contains several charged residues containing an acidic patch on one face of the 
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protein and several basic residues on the other face of the protein. Several of the 
charged residues are located in helices 2, and 3 which correlates well with the 
identified self-association and protein association sites on FAS DD [12] and FADD 
DD [13]. These observations suggest that the mode of interaction of TNFR-1 DD 
5 with TRADD DD may be primarily electrostatic in nature. 

Mutagenesis studies set forth in Example 1 have shown that helix 2 and 
part of helices 3 and 4 contain important residues for both self-association and 
interaction with the TRADD DD. These regions of the TNRF-1 DD protein are also 
correlated with well-defined electrostatic regions on the proteins surface. The 

10 mutants that affect both self-association and interaction with the TRADD DD are 
located on two distinct faces of the protein. One face contains critical residues 
from helix 2 and part of helix 3 (K343A, E344A, R347A, R348A, and D353A). This 
is consistent with the previous observation that the surface formed by a-helices 2 
and 3 has been implicated in the homo- and hetero-association of the death 

15 domains of FAS [12] and FADD [13]. A rotation of 180 0 about the z-axis 
highlights the other surface composed of key residues in helix 4 (E369A and 
Y372F). A similar result was seen with the Tube and Pelle death domains complex 
where (X-helices 4 and 5 of the Pelle death domain interacts with a-helix 6 and the 
unique C-terminal tail of the Tube death domain [19]. 

20 The observation that mutations in two distinct surfaces of TNRF-1 DD 

disrupts both self-association and interaction with the TRADD DD is consistent with 
the proposed mechanism of TNFR-1 DD activity [4]- Since TNRF-1 DD trimer 
formation is a requirement for the recruitment of TRADD DD to the cell surface and 
the interaction of TRADD DD with TNRF-1 DD, the two surfaces identified on the 

25 TNRF-1 DD by site-directed mutagenesis correspond to the self-association binding 
sites (see, Example 1). This would imply that each self-association binding site 
interacts with one TNFR-1 DD in the trimer complex forming a triangular 
arrangement consistent with the X-ray structure of TNFR-1 [49] . Conversely, since 
mutagenesis of either site did not result in the elimination of TRADD DD 
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association with TNRF-1 DD while maintaining the ability of TNFR-1 DD to self- 
associate, neither of the binding surfaces can be exclusively attributed to a TRADD 
DD binding site. Nevertheless, it is still plausible that the TRADD DD binding site 
overlaps with the TNFR-1 DD self-association site. In the case of the FAS [12] and 
5 FADD [13] death domains, mutagenesis data identified distinct binding sites 
associated with hetero-interaction implying a potential difference in the binding 
interaction between TNFR-1 DD and TRADD DD. 

The surface of the R347K TNFR-1 DD structure also contains an acidic 
patch consisting of residues from helices 1 and 6, but mutation of E406A and 
10 E410A (of helix 6) had no effect on self-association or binding to TRADD DD (see, 
Example 1). The assembly of the TNFR-1 DD trimer may present this acidic patch 
as a potential binding site for an undetermined death domain in another signaling 
cascade. Additionally, the acidic patch on TNFR-1 DD may represent the binding 
site of an unidentified regulator of TNFR-1 activity that is potentially released upon 
15 trimer formation. A potential functional role for helix 6 in TNFR-1 DD would be 
consistent with the Tube and Pelle death domains complex where a-helix 6 from 
the Tube death domain comprises part of the interaction site [19]. 

The comparison of TNFR-1 DD with the other members of the death 
domain superfamily of proteins revealed that the major difference in the structures 
20 relative to TNFR-1 DD resides in the lengths and orientation of helices 2, 3, or 4. 
The location of the structural change between the death domains is consistent with 
mutagenesis studies that have identified critical residues for both self-association 
and binding to TRADD DD in these regions of the protein (see. Example 1). 
Particularly, the mutagenesis data identifies the two putative self-association 
25 binding sites for TNFR-1 DD. One site involving 0C2 and 0C3 are consistent with the 
identification of the homo-binding sites for FAS DD [12] and FADD DD [13]. The 
observed differences in the local structure and sequence between TNFR-1 DD and 
the other proteins in the death domain superfamily that correlate with the 
functional regions of the proteins may account for the TNFRl-DD's selectivity and 
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specificity in the TNF signaling cascade while accounting for the similarity in the 
overall fold. The observations that the nature of the protein-protein interaction can 
vary gready further supports the apparendy versatile nature of the death domain 
fold. 
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Table I Structural, Energetic Statistics and Atomic rms Differences 

Structural and Energetic Statistics 



<SA> 



rms deviations from exptl distance restraints (A) (1 167) 
No. of distance restraint violations greater than 0.25A 

rms deviations from exptl dihedral restraints (deg) (92) 
No. of dihedral restraint violations greater than 3 (deg) 

rms deviations from idealized covalent geometry 
bonds (A) 
angles (deg) 
impropers (deg) 



0.0510±0.05 
0 



0.0030±0.0002 

0.409±0.015 

0.3732±0.0123 



(SA)r 



0.0185±0.001 0.0107 
0 0 



0.0655 
0 



0.0015 
0.2860 
0.2477 



Energetics 



E rcMl (kcal mol") 



repel 



E N0E (kcal mol' ) 



E^ (kcal mol") 
E^ (kcal mol") 
Efkcal mol") 



E (kcal mol") 



19±3 
20±3 

0.03±0.05 
13.6±2.5 
15.7±1 
67.1 ±5 



9.6 

6.7 

0.03 

3.2 

6.9 

32 



Non-hydrogen Atomic rms Differences (A) 



All residues 3 
backbone atomsall atoms 



Secondary structure b 
backbone atoms all atoms 



<SA> vs SA 
<SA> vs (SA)r 
(SA)r vs SA 



0.62±0.09 

0.71±0.1 

0.34 



1.50±0.11 
1.70±0.23 
0.81 



0.40±0.07 
0.44±0.08 
0.21 



1.08±0.08 
1.24±0.10 
0.59 



The NMR structures are denoted as follows: <SA> are the final 22 ensemble structures; SA is the 
mean structure obtained from averaging the cartesian coordinates of individual ensemble 
members; and (SA)r is the minimized average structure obtained by regularization of SA. E re , 
was calculated using a final force constant of 4.0 kcal mol" A" 4 with van der Waals hard sphere 
radii scaled by 0.78. E N0E was calculated using a square-well potential with center-averaging and 
a force constant of 50 kcal mol" 1 A'lE^ was calculated using a force constant of 200 kcal mol" 
rad ; E bo^ and E imp were calculated using force constants of 1000 kcal mol"' A' 2 , 500 kcal 
mol rad and 500 kcal mol" rad" 2 , respectively. "In all atomic rms differences calculations, only 
the backbone atoms (N, C a and C) are included in the least squared best fitting. "Core four 
helices. 
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All publications mentioned herein above, whether to issued patents, 

10 pending applications, published articles, or otherwise, are hereby incorporated by 
reference in their entirety. While the foregoing invention has been described in 
some detail for purposes of clarity and understanding, it will be appreciated by one 
skilled in the art from a reading of the disclosure that various changes in form and 
detail can be made without departing from the true scope of the invention in the 
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